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ABSTRACT 

We describe observations at frequencies near 15 GHz of the second 2°x2° field 
imaged with the Cambridge Cosmic Anisotropy Telescope (CAT). After the removal 
of discrete radio sources, structure is detected in the images on characteristic scales 
of about half a degree, corresponding to spherical harmonic multipoles in the range 
i « 330-680. A Bayesian analysis confirms that the signal arises predominantly from 
the cosmic microwave background (CMB) radiation for multipoles in the lower half of 
this range; the average broad-band power in a bin with centroid t — 422 {6 w 51') is 
estimated to be AT/T = 2.ll°;^ x 10"^ For multipoles centred on ^ = 615 [6 « 35'), 
we find contamination from Galactic emission is significant, and constrain the CMB 
contribution to the measured power in this bin to be AT/T < 2.0 x 10~^ [la upper 
limit). These new results are consistent with the first detection made by CAT in a 
completely different area of sky. Together with data from other experiments, this new 
CAT detection adds weight to earlier evidence from CAT for a downturn in the CMB 
power spectrum on scales smaller than 1°. Improved limits on the values of Hq and 17 
are determined using the new CAT data. 
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1 INTRODUCTION 

Observations of spatial fluctuations in the cosmic microwave 
background (CMB) radiation are fundamental to our un- 
derstanding of structure formation in the universe as they 
mark the earliest observable imprints of massive gravita- 
tional structures (see e.g. review by White, Scott & Silk 
1994). The distribution and amplitude of anisotropies in the 
CMB sky over scales from degrees to arcminutes can be 
used to discriminate between competing cosmological the- 
ories. On scales of 0.2°- 2°, inflationary models predict that 
increased power should be seen in the CMB sky due to scat- 
tering of photons during acoustic oscillations of the photon- 
baryon fluid at recombination. Detection and study of these 
acoustic or 'Doppler' peaks in the power spectrum is one 
of the primary goals of CMB astronomy. Furthermore, the 



amplitudes and angular scales of the acoustic peaks pro- 
vide powerful constraints on basic cosmological parameters 
including Hq and Q.. 

The flrst clear indication of a downturn in the power 
spectrum on sub-degree scales was provided by the detec- 
tion of CMB power by the Cambridge Cosmic Anisotropy 
Telescope (CAT) on scales of about half a degree (Scott 
et al. 1996; Paper I). The CAT is a three-element inter- 
ferometer operating at frequencies between 13 and 17 GHz 
(Robson et al. 1993). It is sensitive to structure on angular 
scales of about 10' to 30' over a fleld of view covering 2° x2° 
(primary beam FWHM). This paper describes observations 
of a second field observed with CAT and the detection of 
CMB anisotropies within it, at levels consistent with mea- 
surements in the first fleld. 
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2 OBSERVATIONS AND DATA REDUCTION 

Observations have been made with the CAT of a blank field 
centred at the position 17 00 00 +64 30 00 (B1950), which we 
call 'CAT2'. The field was chosen to be relatively free from 
strong radio sources at frequencies up to 5 GHz (Condon, 
Broderick & Seielstad 1989), lie at high Galactic latitude 
(6 > 30°) and away from known Galactic features (e.g. the 
North Polar spur). 

In periods during the interval 1995 March - 1997 June, 
CAT observed the CAT2 field at three frequencies, 13.5, 
15.5 and 16.5 GHz. The three baselines of the array were 
scaled with frequency to achieve the same resolution at each 
frequency. The resulting synthesised beams measured 20' 
FWHM in right ascension and 24' in declination. The pri- 
mary beam of the telescope, due to the symmetric nearly- 
Gaussian envelope beams of the three horn-reflector anten- 
nas, has a FWHM of 1.96° at 15.5 GHz, scaling inversely 
with frequency. The telescope observes in two orthogonal lin- 
ear polarisations (which rotate on the sky as it tracks), and 
has a system noise temperature of 50 K. An observing band- 
width of 500 MHz was used. Amplitude and phase calibra- 
tions were carried out daily with observations of Gas A (us- 
ing the flux scale of Baars et al. 1977), and cross checked pe- 
riodically by observations of other 15-GHz calibrators from 
the VLA list (Perley 1982). Typical uncertainties in flux scal- 
ing are less than 10%. Observations were generally carried 
out at night (about 80% of the data) or pointing > 90° away 
from the sun to avoid possible solar interference; no extra 
emission from the moon was detected at this declination. 

The CAT2 data were reduced using the same method as 
for the first field, CATl, as described in detail by O 'Sullivan 
et al. (1995). Phase rotation and flux calibration were ap- 
plied first, and the data were then edited and analysed us- 
ing standard tasks in AIPS. Excessively noisy data (reflect- 
ing periods of poor weather) were excised — the visibil- 
ity amplitudes were all checked by eye for periods where 
they regularly exceeded the mean value by more than 3(t 
and then these data (±1 hour) were removed from all base- 
lines. Across the remaining dataset, individual visibilities 
with amplitudes exceeding the 3(j threshold were also ex- 
cluded. In total, about 40% of the data were excluded by 
this process leaving 370, 310 and 1340 hours of good data at 
13.5, 15.5 and 16.5 GHz respectively. Since the atmospheric 
coherence time is very short (10s) compared with the total 
integration time, any remaining atmospheric signals will be 
distributed uniformly across the synthesised map as noise, 
unlike true sky signals which will be modulated by the en- 
velope beam pattern. The efficacy of atmospheric filtering 
for the CAT interferometer has already been demonstrated 
(Robson et al. 1994). CAT's sidelobe response and lack of 
crosstalk and correlator offsets are discussed in O 'Sullivan 
et al. (1995). No radio interference was seen. 



3 SOURCE SUBTRACTION 

Radio sources contributing significantly to the CAT2 image 
were identified and monitored at 15 GHz using the Ryle 
Telescope (RT). Five RT antennas were used in a configura- 
tion giving a synthesised beam of 30" FWHM. The RT has 
an instantaneous field of view of 6' FWHM, but for these ob- 
servations is used in a rastering mode which covers 30' x 30' 



in 12 hours to a typical flux sensitivity of 1.5 mjy/beam. To 
detect sources within the the CAT2 field, the central 2° x 2° 
area was scanned with the RT in raster mode in sixteen 
days. A source list was then compiled, including sources 
listed in the Green Bank 4.85-GHz survey (Condon et al. 
1989). Pointed observations with the RT were then made, 
and repeated regularly to check for variability over the whole 
period of the CAT observations. In all, twenty-nine sources 
were detected, the faintest having a flux density of 4.5 mjy 
at 15 GHz. Flux densities at 13.5 and 16.5 GHz were ex- 
trapolated using spectral information obtained from lower 
frequency surveys (i.e. 4.85 GHz, Condon et al. 1989, and 
1.4 GHz, Condon et al. 1998) where available. A flat spec- 
trum between 13.5 and 16.5 GHz was assumed for three 
sources without other data. After correcting for CAT pri- 
mary beam attenuation, the corresponding flux densities 
were subtracted as point sources from the visibility data. 
Totals of 33mJy at 13.5, 27mJy at 15.5 and 20mJy at 
16.5 GHz were subtracted. The robustness of the source- 
subtraction procedure is illustrated in O'SuUivan et al. The 
strongest source in the field, at position 16 45 32 -|-63 35 29 
(B1950), was highly variable (by factors of up to two over 
periods of a few days) and so was monitored regularly with 
the RT; no residuals at the source position remain after sub- 
traction of the variable source. The successful subtraction of 
this source, l°away from the pointing centre and clearly vis- 
ible at the same position in all three maps, shows that CAT 
phase calibration and pointing are accurate. 



4 THE IMAGES 

The resulting source-subtracted images each show excess 
signal in the central 2° x 2° falling away at larger radii as 
expected from the antenna envelope beam. These are dis- 
played in Figure |l[ The instrumental noise levels (measured 
directly from the visibilities) for each source-subtracted im- 
age are 6.1, 6.5 and 3.5 mJy/beam rms for 13.5, 15.5 and 
16.5 GHz respectively. At the three frequencies, excess pow- 
ers above the intrinsic noise level in the central 2° x 2° area 
of 7.8 ± 1.0, 8.2 ± 1.0 and 8.3 ± 0.5 mJy/beam rms were 
found for the source-subtracted 13.5, 15.5 and 16.5 GHz im- 
ages respectively. Checks were made by splitting the data 
in time (consistent excesses were seen), polarisation (no ex- 
cess power was visible on polarisation difference maps) and 
cross-correlating the source-subtracted image with a recon- 
structed map of the radio sources as in O'SuUivan et al. (no 
correlation was found); most importantly, maps were made 
by correlating orthogonal polarisations and these gave the 
same noise levels as those above. 

To attempt to remove the instrumental response and 
thereby illustrate the distribution of features on the sky, we 
have co-added the data from the three frequencies weighted 
as and CLEANed the final image. This is shown in Figure 
^ and shows the presence of significant features in the central 
region as well as the diminution of the telescope sensitivity 
in accordance with the envelope beam. 

The strongest feature in the CAT2 images is a negative 
one centred at position 17 05 29 -1-64 47 37 (B1950) reach- 
ing —39 mJy at 16.5 GHz (i.e. 5(t, relative to the rms excess 
power in the sky). For comparison, the strongest positive fea- 
ture in the 16.5-GHz dirty map has a peak flux density of 26 
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mjy at position 16 55 55 +63 49 47 (B1950). The negative 
feature can be seen at all three frequencies (most clearly in 
the 13.5 and 16.5-GHz images in Figure both before and 
after source subtraction and in different time cuts, and its 
spectrum is consistent with that of CMB radiation. For ex- 
ample, even before source subtraction the hole reaches —30 
mJy in the 16.5-GHz image, and only three weak sources 
(S16.5 < 10 mJy) lie within 40' of it. No obvious Galactic 
structures at the position of the negative feature were visi- 
ble in IRAS 100 /im (Wheelock et al. 1994) or Hi 21-cm sky 
survey images (Lockman & Dickey 1995) or the lOOpm map 
of Schlegel, Finkbeiner & Davis 1998. Indeed, there is no 
resemblance at all between any of these images of the CAT2 
region and the CAT2 image itself. 

On the scales sampled by CAT, the negative feature is 
unlikely to have been caused by a Sunyaev-Zel'dovich (S-Z) 
effect (Sunyaev & Zel'dovich 1972; Rephaeli 1995) towards a 
single massive cluster. Any cluster which might produce such 
an effect would have to be nearby (subtending a large angle 
on the sky, filling the CAT beam) and/or very massive to 
produce such a strong signature. The strongest S-Z decre- 
ments measured towards nearby clusters at 15 GHz with 
the RT (Grainge et al. 1996; Jones et al. 1993) are about 
—0.5 mJy on arcminute scales. Observing a similar cluster 
at 2 ~ 0.2 (e.g. gas mass 10^*Mq, with a King profile of 
core radius ~ 300 kpc, truncated at about ten core radii) 
with the larger CAT beam of 30' would produce a similar 
decrement amplitude, due to the balance between the efi'ects 
of beam dilution and sampling cluster gas out to a larger 
radius. In order to maximise the S-Z signal in a 30' beam 
with the minimum gas mass, a nearby cluster (z ~ 0.05) 
with IO^^Mq of gas would be required, which is patently 
not observed. A system at any higher redshift would require 
an even higher mass because of beam dilution. 

A central portion of the CAT2 field has been observed 
serendipitously with ROSAT. The direction of the negative 
CMB feature lies close to the edge of the PSPC field (45' 
off axis) and is partially shadowed by the support struc- 
ture of the PSPC detector. An X-ray source (heavily dis- 
torted by the PSPC point spread function) is found nearby 
(at the position 17 05 49 -1-64 45 50, B1950) but its X-ray 
spectrum does not fit a thermal model for any reasonable 
temperature, and is fit much better by a power law, im- 
plying it is not a cluster (but perhaps an AGN or Galactic 
source). Four clusters (Abell 2246 and two others at z — 0.25 
plus one ai z — 0.44) — and an optically luminous QSO 
(HS 1700-1-6416) — have been detected by ROSAT and lie 
within the central 10' of the CAT2 field (Reimers et al. 1997; 
Vikhlinin et al. 1998). All four clusters are fainter than 10** 
ergs~^ at X-ray energies 0.4-2.0 keV and none is apparent 
in the CAT images. 

Finally, we emphasise that the presence of a 5a fea- 
ture in a single CAT field should not be construed as evi- 
dence for non-Gaussianity of the CMB fluctuations: the side- 
lobe structure of the synthesised beam of the three-element 
CAT is significant, and analysis should be carried out in the 
aperture plane — see Section ^. As a check, we have sim- 
ulated CAT images given standard CDM-based realisations 
of CMB structure (using the actual CMB and Galactic mean 
power measured by CAT), and find that features which ap- 
pear as strong as 5a occur in about 10% of cases. We return 
to this point in Section |^. 



5 DETERMINATION OF THE CMB 
COMPONENT 

Due to the limited range of baselines and resulting sparse 
sampling of the u — v plane by CAT, a statistical likelihood 
analysis was employed to estimate the relative contributions 
of CMB and Galactic components given the three-frequency 
CAT data. A Bayesian likelihood method was used, as in 
Paper I, as described by Hobson, Lasenby & Jones (1995). 
This method uses the complex visibility data directly in the 
calculation of the likelihood function to avoid the problem 
of the long-range correlations present between different res- 
olution elements in the image plane. 

As described in Paper I, power was estimated in two in- 
dependent annular bins centred on spherical harmonic mul- 
tipole values of i = 410 and i = 590 and with widths equiv- 
alent to the diameter of the antenna function, thus together 
spanning the range £ « 330 — 680. The noise-weighted cen- 
troid positions for data in each bin are £ — 422 and £ = 615. 
CMB and Galactic signals were modelled as independent 
Gaussian distributions with a power-law spectrum {S oc v"), 
the CMB with a fixed spectral index of +2 in flux density 
and the Galactic spectral index variable between and — 1 
(as expected for Galactic free-free and synchrotron emis- 
sion). 

After marginalising over the Galactic parameters, this 
analysis confirmed that the bulk of the power in the 
16.5 GHz map (Fig. |l|) arises from CMB fluctuations. The 
CMB component was clearly distinguished from possible 
Galactic contamination in the £ = 422 bin; AT/T = 
2.1^0 5 ^ 10~^ was estimated for the CMB signal, compared 
with only AT/T = O.Sto l x 10"^ for any Galactic compo- 
nent. The uncertainties quoted are la values. The CMB- 
Galaxy separation was less certain in the £ = 615 bin — a 
Galactic contribution of AT/T = l.O+i;^ x 10"^ was given 
by the likelihood analysis, which is equivalent to an upper 
limit for the marginahsed CMB power of AT/T < 2.0 x lO"'' 
(la). These values are plotted in Figure ^ The average val- 
ues of AT/T in the CAT2 field agree with the CATl resuh 
(Paper I) within la; for comparison the measured CMB 
powers in the CATl field were AT/T = l.dtol x 10"^ 
{£ = 420) and AT/T = 1.8t°ol x 10"^ = 590). ' 

We have also investigated how much of the CMB power 
is not associated with the strong negative feature discussed 
in Section 4. We removed the dip as a point source from 
the visibilities and repeated the above analysis. Significant 
power remained, at around half of the level with the negative 
feature included. 



6 ESTIMATION OF COSMOLOGICAL 

PARAMETERS USING NEW CAT RESULTS 

The CAT2 points, taken in conjunction with the res ults from 
the Saskatoon experiment (Netterfield et al. 1997), provide 



further evidence for a downturn in the CMB power spec- 
trum for ^^300. To assess the implications for cosmological 
parameters using current CMB data sets, we have extended 
the analysis presented previously by Hancock et al. (1998). 
This analysis used a statistically independent subset of the 
current data and carried out fitting for a range of cos- 
mological models and parameters. The extensions carried 
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Model 




Best Fit 


Marj 


;inalised 


ranges 






value 


/-I 
(1 




(20-) 


Flat models 


h 


0.3 


0.3 - 


0.55 


0.3 - 0.8 




n 


0.88 


0.8 - 


1.15 


0.65 - 1.3 


Lambda models 




0.3 


0.4 - 


- 0.7 


0.3 - 0.7 




h 


0.35 


0.3 - 


0.55 


0.3 - 0.8 




n 


0.9 


0.8 - 


- 1.1 


0.7 - 1.2 


Open models 




0.7 


0.5 - 


- 1.0 


0.3 - 1.0 




h 


0.6 


0.3 - 


- 0.7 


0.3 - 0.8 



Table 1. Results of the fitting analysis. Non-stated parameters 
were marginalised over to estimate the la and 2a ranges. 



and h = 0.6. In conjunction with the CAT points, new re- 
sults forthcoming from the Python, Viper and QMAP^] ex- 
periments may soon be very significant in ruling out open 
models and in further delimiting the Doppler peak, sharp- 
ening up these parameter estimates. Recent OVRO results 
at ^ ~ 590 (Leitch, private communication) agree well with 
the values found by CAT. The joint CAT and OVRO results 
will clearly be very significant in constraining the latest cos- 
mic string power spectrum predictions, which now include 
a cosmological constant (Battye et al. 1998). These predic- 
tions succeed in recovering a significant 'Doppler peak' in 
the power spectrum, but have peak power for £ in the range 
500 to 600, at variance with the trend of the current exper- 
imental results. 



out for the present work were (a) the inclusion of the CAT2 
point at ^ = 422; (b) the inclusion of new points from ex- 
periments Python (Python III, Piatt et al. 1997), MSAM 
(the 2nd and 3rd flights, Cheng et al. 1996, 1997, Ratra et 
al. 199 7), ARGO (Aries+ Taurus regio n, Masi et al. 1996) , 
FIRS ( ICanga et al. 1994| ) and BAM ( [Tucker et al. 1997| ), 
and with the latest calibration correction to the Saskatoon 
data (i.e. increased by 5%, Leitch, private communication); 
(c) consideration of a wider class of cosmological models, all 
treated using exact power spectra rather than the generic 
forms assumed in Hancock et al. (see also Rocha 1997; Rocha 
et al. in preparation). The formalism and approach are oth- 
erwise the same as in Hancock et al. (1998), to which the 
reader is referred. The cosmological models considered were: 

(i) Flat models with A = 0, and with a range of spectral 
tilts (i.e. n 7^ 1); 

(ii) Flat models with A 7^ 0, and a range of spectral tilts; 

(iii) Ope n models with A = and open-bubble inflation 
spectrum (|Hu fc Sugiyama 199*5 ; Ratra & Peebles 1994 
Kamionkowski et al. 1994[ J! 



In cases (i) and (ii ), nucleosynthesis constraints, with 
0.009 < Qth'^ < 0.02 (Copi, Schramm fc Turner 1995) were 



assumed. Theoretical power spectra from Seljak & Zaldar- 
riaga (1996) were used in cases (i) and (ii), and kindly pro- 
vided by N. Sugiyama for case (iii). The parameters fltted 
for (unless flxed) were the Hubble constant h (in units of 
100 kms~^ Mpc~^ ), the spectral tilt n, the cosmological 
constant S^a and the matter density flm in units of the crit- 
ical density. The fiat models are defined by Qm + S7a = 1. 
Note the ranges of h and Qa considered were 0.3 — 0.8 and 
0.3 — 0.7 respectively. 

The results are displayed in Table |^. In each case the 
best fit values of the parameters are shown, together with 
marginalised error ranges. These marginalised errors are 
±1(7 confidence limits formed by integrating over all the 
other parameters with a uniform prior. 

In common with other recent work on fitting to current 
CMB data (e.g. Lineweaver & Barbosa 1998) a tendency to 
low Ho values (except in the case of open models) is found, 
although it is clear from the marginalised ranges that the 
statistical significance of this is not yet high. For compar- 
ison with the results of Webster et al. (1998), who worked 
jointly with recent CMB and IRAS large-scale structure re- 
sults, we note that a flat A model with normalization fixed 
to the COBE results and n = 1 yielded a best fit of = 0.7 



7 CONCLUSIONS 

(1) We have imaged a 2°x2° patch of sky at 13-17 GHz 
with the CAT; this is the second field observed by CAT. 

(2) Signiflcant CMB anisotropy is detected in the CAT2 
field with an average power of AT/T = 2.1 j;";^ x 10"^ for the 
£ = 422 bin, and with an upper limit of AT/T < 2.0 x 10"^ 
for £ — 615 (due to Galactic contamination). 

(3) This new result is consistent with the first detection 
made by CAT in a different area of sky. 

(4) Together with other CMB data over a range of angular 
scales, the inclusion of the new CAT2 detection restricts the 
likely values of cosmological parameters. 
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Figure 1. Images of source-subtracted CAT2 data at a) 13.5, 
b) 15.5 and c) 16.5 GHz. RA and Dec (B1950) and grey scale 
flux ranges (mjy/beam) are indicated. Instrumental noise levels 
are 6.1, 6.5 and 3.5 mJy/beam rms for a), b), c) respectively, 
reflecting the total integration times. 



Figure 2. CLEANed image of co-added CAT2 data (three fre- 
quencies weighted as .) 



Figure 3. Plot of average marginalised power against 
i) in CMB fluctuations (solid lines) and Galactic emission (dotted 
lines) in the CAT2 fleld as estimated by the Bayesian method for 
the two {I = 422 and t = 615) bins. Error bars (Icr) are shown. 
The prediction of a standard CDM model (Om = 1, f^A = 0, 
= 0.05; h = 0.5,n = 1, no tensors) is shown (solid curve) for 
comparison. 
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